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COMPARISON OF LANDSAT-2 AND FIELD SPECTROMETER REFLECTANCE 

! SIGNATURES OF SOUTH TEXAS RANGELAND PLANT COMMUNITIES 1 

1 

1 A. J. Richardson, D. E. Escobar, H. V. Gausman, 

i * 

| and J. H. Everitt 1 

♦ 

I - . - 

; INTERPRETIVE SUMMARY 

i 

We tasted tha accuracy of an atmospheric correction method for 
measuring the reflectance of four prominent south Texas rangeland 
plants, using the Earth Resource Technology (LANDSAT) satellite multi- 
spectral scanner and a ground-based Exotech Model 20 spectroradiometer. 
The atmospheric correction method produced LANDSAT reflectance measure- 
ments of rangeland plants as accurately as the ground-based Exotech 
spectroradiometer. 


1 Contribution from the Soil and Water Conservation Research, Science 
and Education Administration, Agricultural Research, USDA, Weslaco, 
Texas. This study was supported in part by the National Aeronautics 
and Space Administration under Contract No. S-53876-AG. 

Physicist, Biological Technician, Plant Physiologist, and Range 
Conservationist, respectively, SEA, USDA, Weslaco, Texas 78596. 
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INTRODUCTION 

The Earth Resource Technology Satellites (LANDS AT-1 and -2) can 
yield high quality data relevant to the spectral reflectivity of the 
earth’s surface. Since LANDS AT-1 was. launched on August 25, 1972, 
efforts have been made to transform LANDSAT multispectral scanner (MSS), 
digital counts (DC) recorded on computer- compatible tapes (CCT) to 
absolute reflectance valu e s of the earth's surface (Rogers and Peacock, 
1973; Herzog and Sturm, 1975) so as to enhance the use of these data 
for earth resources applications. All of the proposed techniques, how- 
ever, require ground measured solar radiometric data to determine the 
solar and atmospheric parameters that are needed in relating LANDSAT 
count rates to reflectance. However, Ahern et al. (1977) have developed 
a method of using dark targets , such as clear lakes , and atmospheric 
radiative transfer theory (Turner et al., 1971) to estimate the needed 
atmospheric parameters without ground measured solar radiometric data. 

We conducted this study to test Ahern’s method. We compared reflectance 
signatures of four prominent south Texas rangeland plants (Gausmac et al 
1977a and b) obtained by LAND3AT-2 MSS and by the ground-based Exotech 
Model 20 spectroradiometer (Learner et al. * 1973). (Trade names and 
company names are included for the reader^' benefit an.’ do not imply .in 
endorsement or preferential treatment of th'? product by the U.S. r-»r :rr- 
ment of Agriculture.) 
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ATMOSPHERIC RADIATIVE TRANSFER THEORY 
The conversion of LANDS AT digital count data in each band to re- 
flectance (R) at the earth T s surface requires the use of the following 
wavelength dependent atmospheric radiative transfer equation (Turner 
et al., 1971; Rogers and Peacock, 1973; Hulstrom, 1974; Herzog and Sturm, 

* 

1975; Ahern et al. , 1977): 

R = — ~ L P ) * 100, (1) 

ET 

where the atmospheric problem (Fig. 1) for determining R consists of 
evaluating each of the variables defined as follows: 

I 

L - total radiance detected by LANDSAT at the top of the atmosphere 
(mw cm 2 sr ^), 

DC * digital count data recorded on CCTs, 

A,B - LANDSAT radiance calibration coefficients, 

T - vertical atmospheric transmittance of radiant energy from the 
earth’s surface to the LANDSAT MSS, 
t - total optical depth of the atmosphere , ' 

tr - Rayleigh optical depth due to scattering by gaseous molecules, 
tm - Mie optical depth due to scattering by aerosol particulates , 
ta - optical depth due to water absorption, 

_2 

£ - total incident solar irradiance at the earth's surface (mw cm ) j 

I 

(also known as incoming sola r rad iation ; insolation), j 

—2 

Eo - solar irradiance at the top of the atmosphere (mw cm ) 

! (solar constant), 

i 

1 Ts - slant atmospheric transmittance from the sun to the earth's surface 
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-2 

Es - diffuse solar irradiance incident at the earth's surface (raw cm ) , 

_o 

Ed - direct solar irradiance incident at the earth's surface (mw cm ), 

Lp - path radiance detected by LANDS AT at the top of the atmosphere 

, -2 - 1 . ’ . 

C mw cm sr ), 

LI - total radisnce'over a dear lake detected by LANDSAT (mw cm’- sr" 1 ) 

-2 -1 

Lv - radiance from a clear lake water volume Cmw cm sr ), 

-2 -1 

i'Ls - radiance from a clear lake water surface Cmw cm sr ) , 

-2 -I 

Lg - radiance from sun glint due to wave action (mw cm sr ) , 

Rb - background reflectance, 

R - reflectance at the earth's surface, and 
z - solar zenith angle. 

Radiance Detected by LANDSAT CL) 

The first step for solving the atmospheric problem is to convert 

, -2 - 1 , 

the DC data recorded on LANDSAT CCT to radiance (L, mw cm sr ) as 
detected by the LANDSAT MSS at the top of the atmosphere. The equation 

for this operation is as follows: 

✓ 

Li - A^ DC. + B^, where i = LANDSAT band numbers 4, 5, 6, or 7. (2) 

Table 1 lists the LANDSAT radiance calibration constants (A and B) that 
are used for equation (2). 

Atmospheric Transmittance (T) 

The vertical atmospheric transmittance (T) from the L*NDSAT MSS 
sensor to the earth’s surface is computed as: 

T s EXP (-t sec(sensor zenith angle)), (3) 
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where t is the total optical depth of the atmosphere. Even though the 

LANDSAT MSS sensor scans over a range of zenith angles from -5.78 to 

5.78 degrees from the sensor's nadir (Kaneko and Engvall, 1977), it is 

■ 

usually assumed that the sensor zenith angle is zero (vertical). Thus 

for LANDSAT: 

T * EXP (-t). (4) 

Therefore, to be able to calculate T, we only need to know t 
that is a measure of the atmospheric attenuation of incident solar 
irradiance due to scattering and absorption. Scattering effects are 

♦ 

generally assumed to be due to gaseous molecules (tr; Rayleigh optical 
depth) and aerosol particulates (tm; Mie optical depth) (Turner et al., 
1971; Turner and Spencer, 1972). The optical depth due to water absorp- ! 

i 

tion (ta) is assumed to be negligible in LANDSAT bands 4, r . . and 6 but 
not in band 7 (Pitts et al., 1974). Total optical depth (t = tr + tm + 
ta) can be directly measured using a solar radiometer (Rogers and 
Peacock, 1973); however, we used Ahern's et al. (1977) method where t 
is related to Lp through atmospheric radiative transfer theory, using 
a phase function approximation of atmospheric scatterers given by 
Turner et al. (1971) and Turner and Spencer (1972). 

Total Incident Solar Irradiance (E) 

-2 i 

Total incident solar irradiance (E, mw cm ) at the earth’s surface 
may be directly measured with a solar radiometer as Rogers and Peacock 
(1973) and Hulstrom (1974) have shown, or it may be calculated using 
radiative transfer theory (Ahern et al., 1977). 
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As a first step to calculating E, it is necessary to know the solar* 
- 2 ) 

irradiance (Eo, mw cm for each LANDSAT band at the top of the atmos- 

A 

phere such as compiled by Thekaekara at al. (1969) and Thekaekara (1974) 
for the standard earth-sun distance (Table 2). The earth-sun. ratios 
given for each day of the year in ephemeris tables, could be used to 
further refine values of Eo because these values change by 7% annually 
with earth -sun distance. 

Once Eo is known, then the direct incident solar irradiance (Ed, 

_2 

mw cm ) at the earth's surface, as measured with a solar radiometer 
(Rogers and Peacock, 1973), can be computed as: 

Ed * EoTs cos (solar 2enith angle), (5) 

where the slant atmospheric transmittance from the earth's surface to 
the sun (Ts) is 

Ts * EXP (-t sec(solar zenith angle)). (6) 

The solar zenith angle is known for each LANDSAT overpass date. 

_2 

Diffuse incident solar irradiance (Es, mw cm ) at- the earth's 
surface, also known as skylight, may be measured by shadowing a solar 
radiometer detector. For this study, Es was calculated using the phase 
functions of atmospheric scatterers as given by Turner et al. (1971). 

Therefore, once Ed and Es are known then the total incident solar 
irradiance at the earth's surface is calculated as: 

E s Ed + Es. (7) 

Path Radiance (Lp) 

-2 -1 

Path radiance (Lp, mw cm sr ) is difficult to determine because 

it cannot be rsasured directly. It depends on a complex interaction 
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between atmospheric scattering and absorption of incident solar irra- 
diance and reflected solar radiance from background albedo (Turner, 1975) 

that is scattered into the optical path of the LANDSAT MSS. Thus, 

* «• 

several methods have been proposed to-.inferpath radiance indirectly. 

Ground-based- solar radiometric measurements of diffuse sky irra- 
diance have been used to indirectly derive path radiance using methods 
given by Gordon et al. (1973), Rogers and Peacock (1973), and O'Neill 
and Miller (1977). In addition, Hulstrom (1974) used a plot of L 
against ground-based measurements of reflectance for various naturally 
occurring calibration targets on the earth : s surface to determine oath 
radiance. Such a plot does not pass through the origin; instead at 
zero reflectance, Lp * L. The weakness of these methods is that they 
depend on ground-based solar radiometric measurements that are not 
readily available. 

Ahern et al. (1977) used the radiance of dark targets, such as 
clear lakes, to determine path radiance. He found that the radiance 
aver a clear lake (LI), at the top of the atmosphere, is the sum of 
several terms : 

LI = (Lv + Ls + Lg) T + Lp (8) 

•fhere Lv is the radiance from the water volume, Ls is the radiance from 
the water surface , and Lg is the radiance from sun glint due to wave 
iction caused by high winds or solar zenith angles less than 30°. From 
lata given by Ahern et al. (1977), Lv = RvE, Ls » 0.00$ Es, and Lg =-0. 
U.so, Rv was estimated from Ahern's data using the following empirical 
squat ion : 
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Rv « 0.0035 - 0.0036X, (9) ; 

where the wavelength (X) ranges from 0.4- to 3.0-yra. Thus, path radiance 
was calculated as: ; 

*. # • i 

Lp « LI - RvET - O.OOSEsT. (10) ' 

I 

Reflectance Variation With Sun Angle 

The LANDSAT MSS and ground-based spectroradiometer measured reflec- 
tance at the earth's surface at different solar zenith angles for the 
same plant. Smith et al. (1975) and Duggin (1977), found that LANDSAT 
reflectance signatures may need to be corrected for plant canopy reflec- 
tance variations with sun angle. However, Lemme and West in (1978) 
observed that reflectance data collected from about 1015- to 1500 -h CDT 
show minimal effect due to sun angle variation. As a result, we did not 
attempt to make any sun angle corrections, because the data for both 
LANDSAT and the ground-based spectroradiometer were collected within 
this time range. 

Background Reflectance (Rb) 

An estimate of the average background reflectance (Rb) is needed 
when using the nhasu function anproximations of atmospheric scatterers 
that relate Lp to t. We used Ahern's et -oil. (1977) approach, which 
calculates Rb with the following equation: 

Rb 3 — (11) 

Eo cos (solar zenith angle) 
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The value for L vat determined by averaging the LANDS AT DC values from 
a 512 by 512 pixel matrix for a study area of interest and then using 
the A and B values in Table 1 to convert to mean radiance (L). 

EXPERIMENTAL PROCEDURES 

Four prominent rangeland plant communities in south Texas are 
(Kuchler, 1964; Davis and Spicer, 1965): Ci) live oak C Cuercus 

virginiana Mill.), a tree that grows on deep sands in formations ranging 
from dense, uniform stands to frequent thickets or motts in underbrush; 
Cii) silver leaf sunflower ( Helianthus argophyllus Torn, and Gray), a 
taprooted annual weed that has white-tomentose plant parts, germinates 
in April or May, reaches leaf pubescence peak in July, and flowers in 
late summer or fall; (iii) cenizo (Lcucophyllum frutescens (JBerland) 

I. M. Johnst.), a woody shrub that grows as either dense or sparse 
stands among a wide variety of woody shrubs on shallow soils; and Civ) 
hoiey mesquite ( Prosopis glandulosa Torr.) that grows as motts or dense 
stan’s on a variety of soil types Cdeep sands, sandy loams, clay loams, 
or heavy clays). 

We used LANDSAT MSS CCT and corresponding color images Cl; 1,000, 000 
scale) for a LANDSAT-2 overpass on June 2, 1977 (Scene I.D. 2862*16000) . 
All four of the LANDSAT KSS bands were used, covering the 0.5- to 1.1-ym 
spectral region . This overpass provided DC data for a 185- by 185-km 
scene that included sample sites, near Sarita, Alice, and Edinburg, Texas 
for the four plant communities. 
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We averaged LANDSAT MSS DC data over 417 training pixels (picture | 

I 

elements) collected from the four plant community sample sites and a 

clear lake. The average of the DC values from the clear lake (ID was 

• # 

used to estimate Lp from equations ‘<8) , (9), and (10). The average of 

the pixels within a 512 by 512 pixel area near Sarita, Ttxas, was used 

to estimate KET from equation (11). Then the Lp and RET averages were 

used with Ahern's method to calculate the T and E which were used with 

Lp In equation (1) to convert the LANDSAT-2 DC averages for each plant 

community to plant reflectance at the earth's surface. 

The field reflectance snectra were previously collected by Gausman 

at al. (1977a) for the sllverleaf sunflower and by Gausman at al. (1977b) 

for the live oak, cenizo, and honay mesquite, over the 0.5- to 2.5-ym 

waveband, during tha 1976 growing season with a Exotech Model 20 spectro- 

radiometer (Learner et al., 1973). The sensor had a 15 -degree fie^d-of- 
2 

view (C.5 m ) and was placed 3- to 3.4-m above each of fiva randomly 
selected canopies for each plant community sample site. 

Using correlation techniques, we analyzed the reflectance data from 
both LANDSAT and spec troradiometer sensors, at the mid-band wavelength 
intervals of the LANDSAT MSS (0.55-, 0.65-, 0.75-, and 0.95-pra). Such a 
correlation will have unit slope and zero interceot if the measured 
reflectance from both sensors for the same plants were identical. There- 
fore, we conducted a t-test analysis to test for a significant deviation 
of the slope from unity end of the interceot from zero. 
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*JCFERIKENTAL RESULTS AND DISCUSSION 
The value of radiance over the clear lake (LI) for band 4 (Table 2) 
waa high (0.461 mw cm * *r A ) as conpared with Ahern's et al. (1977) 

— T * • * 

average value of 0.329 mw cm * sr . This high value overestimated the 

“2 *1 

path radiance (0.438 raw cm sr ) that was used to determine atmospheric 
optical depth (t). Thus, for band 4, t * 0.791, which corresponds to 
horizontal visible range of only 10 km (Potter and Shelton, 1974). The 
horizontal visible range on June 2, 1977 near the rangeland sites was 
probably more than 23 km. Probably the lake we used as a clear water 
reflectance standard was more turbid than we originally assumed. Also, 
the Turner model probably calculates too little path radiance for a 
given optical depth. The Ip value for band 5 was not overestimated as 
much as that for band 4. Values for bands 6 and 7 seemed reasonable as 
compared with Ahern's data. 

The solar and atmospheric parameters given in Table 2 were used 
to convert the LANDSAT-2 digital count data in Table 3 to plant reflec- 
tances for the four rangeland plant communities using equation (1). 

The four plant communities were ranked In descending order by their 
reflectance values in LANDSAT band 7 so that values for sunflower > 
live oak > mesquitc > cenizo. This ranking agreed with previous reflec- 
tance results using ground-based spectroradiometer measurements collected 
by Gausman et al. (1977b) for the three woody canopies, but it differed 
from their reported leaf ground cover values, where the ranking was live 
oak > cenizo > mesquite. The silverleaf sunflower's white -tomentose 
condition apparently caused its reflectance to be higher than that of the 
woody plants. H 
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Figure 2 compares the LANDSAT reflectance values (.) from Table 3 
with tne previously determined ground-based spectroradiometer reflectance 

measurements (solid lines) for the. same plant communities (Gausman et al. 

* 

1977a,b). The values seem quite compatible, except that the LANDSAT 
reflectance values in bands 4 and 5 for Cenizo and bands 6 and 7 for live 
oak were lower than the corresponding ground-based reflectance measure- 
ments. Apparently, the undetermined amount of live oak vegetation cover 
was not very high so that the reflectance in band 6 and 7 was decreased 
due to integrating more soil end shadow reflectances in with plant 
reflectance over a wide ground area (Richardson et al. , 1975). 

Figure 3 shows the slope and intercept results of correlating the 
LANDSAT and ground-based spectroradiometer reflectance measurements. The 

correlation of the reflectance values between the two sensors was highly 
2 

significant (r = 0.924) and a t-test analysis showed that the slope 
(0.994) did not differ significantly from unity nor did the intercept 
(1.55) differ significantly from zero. Thus, these results indicated 
that the LANDSAT MSS could be calibrated for solar and atmospheric 
variations to yield reflectance measurements at the earth's surface that 
were not significantly different from ground-based spectroradiometric 
reflectance measurements , even though the lake used as a clear water 
reflectance standard may have been somewhat turbid. 
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Ta^la 1. LANDSAT-1 and -2 calibration constants for convincing digital count ratas 
to radianco as measured by tha LAHDSAT multispectrai scannar (MSS) at tha 
top of tha aarth'a atooaphara. (From Pottar (1972), Rogers and Peacock 
(1973), Herzog and Sturm (1975), Ottermaa and Fraaar (1976), and LANDSAT 
Newsletter #15.) 


LAHDSAT 

HSS 

Sensor 

LAHDSAT Calibrat. 

n Constants* 


Life-Span of 
LANDSAT Calibration 
Constants 

MSS4 

MSS5 

MSS6 

■" MS$7 

1 A 

0.0195 

0.0157 

0.0138 

0..0730 

8/25/72 to 1/10/78 

B 

0.0 

0.0 


0.0 


2 A 

0.157 

0.0117 

0.0105 

0.0837 

1/22/75 to 7/15/75 

B 

0.10 

0.07 

0.07 

0.14 


2 A 

0.0201 

0.0134 

0.0115 

0^0603 

7/16/75 to present 

4 1 

B:8foos 

m, 

8:8!l3 


7/31/78 to present 

B 

0.04 

0.03 

0.03 

0.03 



*• Tha radianca units for A and B ara mw car” 2 sr”* count”* and nw an” 2 sr”* # 
respectively. 


Table 2. Solar and atmospheric variables determined for a June 2, 1977 LAHDSAT overpass 

(scene I. D. 2862-16000) of rangeland comunities located in south Texas. Solar 
zenith angle was 34 degrees. 


1AITDSAT 

HSS 

Bands 

Clear 

Lake 

Radiance 

(LI) 

Path 

Radiance 

(Lp) 

Solar 

Constant 

(Bo)* 

Diffuse 

Radiance 

(Es) 

Direct 

Radiance 

(Ed) 

Atmospheric 
Transmittance 
And Optical 
Depth 

(T) (t) 

Background 

Reflectance 

CBS) 


-2 -1 
mw cm sr 

-2 -1 
mw cm sr 

-2 

mw cm 

-2 

mw cm 

mw cm” 2 




4 

0.461 

0.438 

17.3 

7.7 

5.5 

0.453 

0.791 

0.133 

5 

0.274 

0.253 

15.1 

5.8 

6.2 

0.5S4 

0.590 

0.118 

6 

0.163 

0.148 

12.4 

2.8 

7.3 

0.751 

0.285 

0.213 

7 

0.170 

0.155 

25.1 

3.0 

17.7 

0.872 

0.136 

0.242 


* From Thekaekara at al. (1969), Rogers and Peacock (1973), and Otterman and Fraser (1976). 

Table 3. Digital count (DC) data, radiance at top of atmosphere (L), and reflectance 
(R) measured by LAHDSAT- 2 on June 2, 1977 (scene I. D. 2862-16000) for four 
typical rangeland vegetation communities. Solar zenith angle was 34 degrees. 


Rangeland LAHDSAT MSS Bands 

Vegetation — - — — — — - ——————— — ■ 


Communities 

DC4 

DCS 

DC6 

DC7 

L4 

L5 

L6 

L7 

R4 

RS 

R6 

R7 






- - mw cm” 2 

sr” 1 


% 

% 

% 

% 

Silverleaf 

Sunflower 

26.2 

25.6 

67.9 

32.8 

0.61 

0.40 

0.84 

2.09 

8.8 

7.1 

28.9 

33.7 

Live Oak 

22.0 

21.3 

$6.0 

23.3 

0.52 

0.35 

0.70 

1.82 

4.4 

4.4 

23.2 

29.0 

Mesquite 

23.8 

25.5 

51.6 

24.7 

0.56 

0.40 

0.65 

1.60 

6.3 

7.1 

21.1 

25.2 

Cenizo 

21.1 

20.2 

47.9 

23.1 

o.so 

0.33 

0.61 

1.50 

3.5 

3.7 

19.3 

23.5 


pane IS 



Figure 1 


Figure 2 


Figure 3 


FIGURE LEGENDS 

t 

Generalized diagram of the atmospheric problem for converting 
LANDS AT digital count data to reflectance (R). The solar and 
atmospheric variables involved were defined previously. 
Comparison of ground-based Exotech spectroradiomatric (solid 
line) and LANDS AT (.) reflectance measurements of four south 
Texas rangeland plants. 

Correlation of ground-based spectroradiomatric (EXOTECH) and 
LANDSAT-2 MSS reflectance measurements at wavelengths 0 . 55- , 
0.65- , 0.75- , and 9 .95 -pm for four south Texas rangeland 
plants . 
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TOP OF THE ATMOSPHERE 






SILVKRLKAP 















R - 1.55 + 0.994 R 
r - 0.961 

Intercept - 1.55 - 3.15 
Slope - 0.994 - 0.076 


Fig. 3 



